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Abstract 

Experimental  data  were  collected  for  the  purpose  of  validating  multiphase 
combustion  models  and  submodels.  A spray  combustor  was  fabricated  that  permits 
well-defined  boundary  and  input  conditions,  and  measurements  were  carried  out  that 
characterize  the  fuel  spray,  wall  temperatures,  gas  temperatures  and  species 
concentrations  at  the  reactor  exhaust.  The  operating  conditions  for  the  baseline  case 
were  defined  by  NIST  personnel  and  industrial  collaborators  involved  in  model 
development.  A methanol  spray,  generated  with  a pressure-jet  nozzle,  was  studied  with 
phase  Doppler  interferometry  to  measure  the  size,  velocity,  number  density,  and  mass 
flux  of  the  fuel  droplets.  Fourier  transform  infrared  spectroscopy  was  used  to  measure 
the  species  concentrations  in  the  reactor  emissions.  The  conversion  of  methanol  in  the 
reactor  was  found  to  be  ~ 80  %.  Inlet  and  outlet  air  velocity  (including  inlet  air  swirl 
intensity)  and  heat  flux  measurements  at  the  reactor  wall  are  planned.  Flow  field 
velocity  and  temperature  measurements  are  also  planned  to  complete  the  baseline 
database. 

Key  words:  CFD  model  validation,  multiphase  combustion,  spray  combustion,  phase 
Doppler  interferometry,  FTIR  spectroscopy 

Certain  commercial  equipment,  materials,  or  software  are  identified  in  this  publication  to  specify 
adequately  the  experimental  procedure.  Such  identification  does  not  imply  recommendation  or 
endorsement  by  the  National  Institute  of  Standards  and  Technology,  nor  does  it  imply  that  the  materials 
or  equipment  are  necessarily  the  best  available  for  this  purpose. 
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1.  Introduction 

Control  of  process  efficiency  and  the  formation  of  species  byproducts  from 
industrial  thermal  oxidation  systems  (e.g.,  power  generation  and  treatment  of  liquid 
chemical  wastes),  is  generally  based  on  a priori  knowledge  of  the  input  stream  physical 
and  chemical  properties,  desired  stoichiometric  conditions,  and  monitoring  of  a few 
major  chemical  species  in  the  exhaust.  Optimization  of  the  performance  of  these 
systems  is  relying  increasingly  on  computational  models  and  simulations  that  help 
provide  relevant  process  information  in  a cost-effective  manner.  In  general,  there  is  a 
dearth  of  reliable  data  for  specifying  model  initial/boundary  conditions,  and  a need  for 
validation  of  the  numerical  codes.  System  performance  is  dependent  on  the  liquid 
atomization,  chemistry,  aerodynamic  design,  and  the  degree  of  liquid/air  mixing  within 
the  reactor.  Provision  of  in-situ,  real-time  data  on  the  characteristics  of  the  droplet  field 
and  flame  structure,  and  an  understanding  of  its  interrelationship  with  the  system 
operating  conditions,  heat  transfer,  and  particulate  and  gaseous  byproducts,  is  crucial 
for  the  development  and  validation  of  advanced  computational  models,  diagnostics,  and 
instrumentation. 

Although  computational  fluid  dynamics  (CFD)  offers  a cost-effective  alternative  to 
experiments,  the  accuracy  of  the  CFD  model  must  first  be  assured.  This  should  be 
accomplished  in  two  ways:  verification  and  validation  (Oberkampf  et  al.,  1997). 
Verification  involves  ensuring  that  the  algebraic  and  differential  equations  within  the 
model  have  been  accurately  solved.  In  addition  to  verifying  that  the  numerical  code 
arrives  at  the  correct  solution,  it  is  also  necessary  to  determine  if  the  correct  model  has 
been  solved.  This  is  the  validation  step.  The  objective  of  this  report  is  to  provide 
experimental  data  for  the  purpose  of  CFD  model  and  submodel  validation,  and 
interested  parties  are  encouraged  to  use  the  enclosed  benchmark  spray  combustion  data 
set.  A file  (mesh.txt)  in  PLOTS  D format*  is  provided  in  the  accompanying  diskettes 
that  may  be  used  to  assist  in  initiation  of  the  simulation  of  the  NIST  Reference  Spray 
Combustion  Facility. 

This  report  presents  data  obtained  fi-om  a baseline  spray  flame  within  the  reference 
spray  combustion  facility  at  NIST.  Methanol  was  chosen  as  the  fuel  for  the  baseline 
case  because  the  thermodynamic  and  kinetic  data  necessary  to  model  the  gas  phase 

* The  file  contains  a sample  2-D  structured  mesh  of  the  numerical  simulation  of  the  NIST  facility.  The 
first  row  in  file  is  the  number  of  domains,  and  the  second  row  is  the  number  of  cells  in  x,  y,  and  z 
dimensional  planes  for  each  domain.  The  actual  mesh  coordinates  for  x,  y,  and  z follow  in  sequential 
order  across  the  four  columns.  The  single  digit  numbers  at  the  end  of  the  file  detail  whether  or  not  the 
mesh  node  is  blocked.  (Courtesy  of  CFD  Research  Corporation,  Huntsville,  AL.) 
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combustion  are  readily  available  (Afeefy  et  ai,  1998;  also  available  on  the  World  Wide 
Web  at  http://webbook.nist.gov/).  The  spray  data  presented  Avere  collected  non- 
intrusively  using  phase  Doppler  interferometry  (PDI).  Data  are  presented  for  the  size 
and  velocity  distributions  of  the  fuel  droplets,  the  droplet  number  density,  and  the 
volume  flux  of  fiiel  droplets  within  the  spray.  The  enclosed  combustion  chamber 
provides  well-characterized  boundary  conditions,  and  wall  temperature  data  are 
provided  as  a function  of  axial  position.  Temperature  and  species  measurements 
obtained  at  the  reactor  exit  are  provided,  which  can  be  used  for  boundary  conditions  or 
validation  of  computational  results.  Gas-phase  velocity  measurements  are  planned 
using  PDI  by  seeding  the  combustion  air  with  spherical  salt  particles  with  diameters  of 
order  one  micrometer,  and  these  results  will  be  compared  to  measurements  obtained 
with  a pitot  tube.  Heat  flux  measurements  at  several  elevations  along  the  reactor  wall 
are  also  planned. 

2.  Description  of  Experimental  Facility  and  Procedure 

2.1.  Spray  Combustion  Reactor  - Baseline  Case 

The  experiments  were  conducted  in  a stainless  steel,  enclosed  spray  combustion 
facility,  shown  in  Fig.  lA.  The  experimental  facility  includes  a swirl  burner  v^th  a 
movable  12-vane  swirl  cascade.  The  cascade  is  adjusted  to  impart  the  desired  degree  of 
swirl  intensity  to  the  combustion  air  stream  that  passes  through  a 0.10  m diameter 
passage  and  coflows  around  the  fuel  nozzle.  The  swirl  intensity  is  a measure  of  the 
angular  momentum  of  the  combustion  air.  It  is  characterized  by  the  swirl  number,  S, 
defined  as  the  ratio  of  the  axial  flux  of  angular  momentum  to  the  axial  flux  of  linear 
momentum  (Gupta  et  al.,  1984).  The  vane  angle  and  combustion  air  flow  rate  were 
50°  ± 1°  and  56.7  ± 1.7  m^  h’^  respectively.  Using  the  correlation  suggested  by  Beer 
and  Chigier  (1972),  the  experimental  conditions  used  here  should  result  in  a swirl 
number  of  S = 0.48,  which  is  a swirling  flow  of  moderate  intensity.  However,  recent 
CFD  modeling  of  the  vane  cascade  and  combustion  air  passage  yields  a swirl  number 
of  S = 0.54.  Because  flame  stability  and  the  combustion  process  are  strongly  dependent 
on  the  amount  of  swirl  present,  experimental  measurements  are  in  progress  to 
investigate  the  discrepancy  between  the  above  estimates. 

Figure  IB  presents  a close-up  view  of  the  burner  and  nozzle.  The  liquid  fuel  is 
forced  through  a pressure-jet  nozzle  and  forms  a hollow-cone  spray  with  a nominal  60° 
full  cone  angle.  A parallel  program  is  underway  to  develop  a reference  atomizer  with 
well-defined  droplet  size  and  velocity  distributions.  The  atomizer  will  then  provide 
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known  spray  characteristics  at  the  inlet,  and  define  this  boundary  condition  for  CFD 
models.  This  atomizer  is  planned  to  be  available  in  the  next  phase  of  this  program. 
Methanol  was  used  for  these  experiments,  and  the  flow  rate  was  maintained  at  3.0  ± 
0.02  kg  h‘^  The  fuel  and  combustion  air  were  introduced  into  the  reactor  at  room 
temperature.  The  fuel  flow  rate,  combustion  air  flow  rate,  wall  temperatures,  and 
exiting  gas  temperatures  were  monitored  and  stored  using  Lab  VIEW  data  acquisition 
software  running  on  a 166  MHz  personal  computer.  The  operating  conditions  for  the 
baseline  case  are  summarized  in  Table  1. 

The  burner  is  enclosed  within  a stainless  steel  chamber  to  provide  improved 
reproducibility  and  control  of  the  spray  flame.  The  chamber  height  is  1.2  m and  the 
inner  diameter  is  0.8  m.  Several  windows  provide  optical  access,  and  a stepper-motor- 
driven  traversing  system  translates  the  entire  bumer/chamber  assembly  permitting 
measurements  of  spray  properties  at  selected  locations  downstream  of  the  nozzle. 
Additional  details  on  the  design  of  the  burner  are  available  in  the  literature  (Presser  et 
al.,  1993).  The  relevant  dimensions  necessary  for  modeling  the  facility  are  presented  in 
Fig.  2.  Note  that  the  reactor  exit  is  off-axis,  which  makes  the  problem  non- 
axisymmetric. 

2. 2.  Phase  Doppler  Interferometer 

Since  its  introduction  (Bachalo  and  Houser,  1984a;  Bachalo  and  Houser,  1984b), 
phase  Doppler  interferometry  (PDI)  has  been  used  to  characterize  sprays  in  areas  such 
as  liquid  fuel  spray  combustion,  coal  slurry  combustion,  coatings,  pesticides,  fire 
suppression,  and  others.  PDI  is  an  extension  of  laser  Doppler  velocimetry  that 
measures  droplet  size  as  well  as  velocity.  Phase  Doppler  techniques  involve  creating 
an  interference  pattern  in  the  region  where  two  laser  beams  intersect  and  results  in  a 
region  consisting  of  alternating  light  and  dark  Singes.  The  region  where  the  laser 
beams  intersect  is  called  the  probe  volume  or  sample  volume.  Due  to  the  interference 
pattern,  a droplet  passing  through  the  probe  volume  scatters  light  that  exhibits  an 
angular  intensity  distribution  which  is  characteristic  of  the  size,  refractive  index,  and 
velocity  of  the  droplet.  For  a droplet  with  known  refractive  index,  the  size  and  velocity 
can  be  determined  by  analyzing  the  scattered  light  collected  with  several 
photomultiplier  tubes. 

The  PDI  is  a single-point  (or  spatially  resolved)  diagnostic  instrument  in  that  it 
obtains  information  about  the  spray  at  a single  point  in  space.  Only  by  moving  the 
probe  volume  (or,  equivalently,  the  spray)  can  one  map  out  the  spatial  distribution  of 
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the  spray  characteristics.  The  PDI  is  also  a single  particle  instrument  in  that 
information  is  obtained  for  only  one  droplet  at  a time.  This  offers  advantages  over 
integrating  techniques  because  the  characteristics  of  a particular  droplet  (size,  velocity, 
etc.)  can  be  recorded  and  the  data  can  be  separated  into  classes  (size  classes,  velocity 
classes)  to  further  characterize  the  spray  system. 

The  spray  measurements  reported  here  were  made  using  a two-component  Phase 
Doppler  Particle  Analyzer  from  Aerometrics  to  measure  the  droplet  size  and  velocity 
distributions  and  the  spray  intensity.  The  instrument  is  composed  of  the  following 
components:  (i)  transmitter  (model  XMT-1100-4S),  (ii)  motor  controller  (model  MCB- 
7100-1),  (iii)  receiver  (model  RCV-2100),  (iv)  counter-type  signal  processor  (model 
PDP-3100),  and  (v)  data  management  system  (model  DMS-4000-5).  The  PDI  was 
controlled  using  Aerometrics  PDPA  software  version  4.27J  run  on  a 66  MHz  personal 
computer.  The  receiving  optics  were  aligned  at  a 30°  scattering  angle  measured  from 
the  direction  of  propagation  of  the  laser  beams,  and  the  transmitting  and  receiving 
optics  were  aligned  at  the  same  elevation.  A 5 W Spectra  Physics  argon  ion  laser 
operating  in  multi-line  mode  was  used  as  the  illumination  source.  The  blue 
(wavelength  = 488  nm)  and  green  (wavelength  = 514.5  nm)  lines  of  the  argon  ion  laser 
were  separated  and  focussed  by  the  transmitting  optics  to  intersect  and  form  the  probe 
volume.  Additional  details  of  the  optical  arrangement  are  available  elsewhere  (Presser 
etal,  1994). 

The  location  of  the  reaction  chamber  is  controlled  by  a series  of  stepper  motors  that 
permit  the  reactor  to  be  moved  in  three  directions.  For  the  experiments  described  in 
this  report,  the  position  of  the  reactor  was  varied  in  only  two  directions,  which  enabled 
the  spray  to  be  probed  in  the  axial  and  radial  directions.  The  reactor  was  not  moved  in 
the  direction  of  propagation  of  the  PDI  laser  beams.  The  third  dimension  is  usually 
used  for  alignment  of  the  optics,  but  could  be  used  to  determine  the  tangential 
component  of  velocity  if  necessary.  Table  2 presents  the  uncertainties  associated  with 
determining  the  location  of  the  PDI  probe  volume  (i.e.  the  location  where  the  beams 
intersect)  relative  to  the  spray  nozzle.  The  uncertainty  in  the  axial  (Az)  position,  which 
is  a quantitative  estimate  of  the  maximum  that  this  measurement  could  be  in  error,  is 
primarily  due  to  the  difficulty  in  aligning  the  laser  beams  with  the  surface  (face)  of  the 
nozzle  (z  = 0).  Uncertainties  in  the  two  horizontal  positions  (Ax  and  Ay)  result  from 
the  difficulty  in  locating  the  probe  volume  exactly  above  the  nozzle.  The  uncertainty  in 
the  spray  angle,  A0,  results  from  an  estimate  of  the  maximum  that  the  bumer/nozzle 
assembly  deviates  from  the  vertical  position.  This  uncertainty,  A0,  also  depends  on  if 
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the  orifice  is  designed  and  fabricated  properly  so  that  the  spray  exits  vertically.  This 
deviation  results  in  a negligible  error  in  the  measurements  close  to  the  nozzle,  but  can 
become  significant  further  downstream  as  shown  in  the  expression  for  Ar  in  Table  2. 

To  characterize  the  droplets  as  a function  of  location  within  the  spray,  it  is  first 
necessary  to  establish  the  origin  of  the  coordinate  system.  This  was  done  using  the 
following  procedure:  (i)  align  the  PDI  receiving  optics  with  the  probe  volume  located  in 
a downstream  region  of  the  spray  where  the  number  density  of  droplets  is  high  and  the 
light  scattering  is  intense,  (ii)  align  the  surface  of  the  nozzle  with  the  PDI  probe  volume 
(z  coordinate),  (iii)  adjust  the  horizontal  position  perpendicular  to  the  direction  of 
propagation  of  the  PDI  laser  beams  (x  coordinate),  and  (iv)  adjust  the  horizontal 
position  in  the  direction  of  propagation  of  the  laser  beams  (y  coordinate)  so  that  the  PDI 
probe  volume  corresponds  with  the  center  of  the  spray  nozzle.  This  alignment 
procedure  is  carried  out  visually. 

The  alignment  of  the  PDI  receiving  optics  was  performed  as  described  in  the 
operating  manual,  and  the  details  will  not  be  repeated  (Aerometrics  Inc.,  1987).  After 
the  axial  coordinate  was  adjusted  (step  (ii)  in  the  alignment  procedure),  the  burner  was 
lowered  slightly  (~1  mm)  so  that  the  probe  volume  was  located  just  downstream  of  the 
nozzle.  The  burner  was  then  translated  in  the  x-direction  until  the  probe  volume  was 
directly  above  the  nozzle  face.  Adjustment  in  this  direction  (x  coordinate)  was  more 
straightforward  than  in  the  other  horizontal  direction  (y  coordinate)  because  one  can 
look  down  vertically  at  the  crossing  laser  beams  to  facilitate  the  alignment.  In  the  y- 
direction,  the  probe  volume  is  elongated  due  to  the  small  intersection  angle  of  the  laser 
beams,  and  this  increased  the  difficulty  in  aligning  the  probe  volume  in  this  direction. 
Therefore,  the  alignment  in  the  y-direction  was  accomplished  utilizing  the  viewing  port 
in  the  PDI  receiving  optics,  as  shown  in  Fig.  3.  Figure  3A  depicts  the  image  that 
appears  in  the  viewing  port  when  the  probe  volume  is  located  1 mm  above  the  nozzle  (z 
= 1 mm).  The  horizontal  lines  correspond  to  the  intersecting  laser  beams,  and  the 
vertical  lines  represent  the  100  |im  slit  that  limits  the  amount  of  light  detected  by  the 
photomultiplier  tubes.  Because  of  the  very  high  intensity  of  the  elastic  scattering  at  this 
location,  the  spray  above  the  probe  volume  is  illuminated  sufficiently  to  be  observed 
through  the  viewing  port.  The  cone  formed  by  the  spray  can  be  used  to  align  the  probe 
volume  in  the  y-direction.  If  the  burner  is  lowered  so  that  the  probe  volume  is  ftirther 
downstream,  the  image  appears  like  that  shown  in  Fig.  3B,  which  corresponds  to  z = 3 
mm. 
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2.3.  Fourier  Transform  Infrared  Spectrometer 

Gas-phase  species  concentrations  were  measured  using  Fourier  transform  infrared 
(FTIR)  spectroscopy.  An  FTIR  spectrometer  (BioRad  FTS-40)  equipped  with  a DIGS 
detector  was  used  for  extractive  sampling  of  chemical  species  in  the  combustor 
emissions.  A gas  sampling  system,  consisting  of  an  air-cooled  sampling  probe,  a 
heated  gas  line,  and  a vacuum  pump,  facilitated  the  transport  of  the  sample  gas 
extracted  from  the  spray  combustor  into  the  single  pass  gas  cell  in  a continuous  manner. 
The  sampling  probe  was  designed  to  aerodynamically  quench  chemical  reactions 
occurring  within  the  gasses  being  sampled.  The  sampling  gas  line  was  also  provided 
with  a means  for  purging. 

In  this  investigation,  the  FTIR  spectrometer  was  operated  with  2 cm'^  resolution. 
The  gas  cell,  equipped  with  KBr  windows,  was  unheated  and  had  a 0.1  m path  length. 
The  sample  gas  enters  the  cell  at  one  end  and  continuously  exits  at  the  other  end.  The 
pressure  inside  the  cell  was  kept  steady  at  a constant  sub-atmospheric  level.  During 
operation,  the  sampling  gas  line  was  purged  frequently  with  a stream  of  purified  air. 
The  data  acquisition  and  analysis  functions  were  accomplished  with  the  software 
provided  by  the  instrument  vendor.  In  each  test  run,  care  was  taken  to  allow  sufficient 
time  for  equilibrating  the  flow  and  thereby  assuring  a steady  state  flow  condition. 
Generally,  three  to  five  spectra  were  obtained  at  each  probe  location.  Once  the  data 
acquisition  at  one  location  was  completed  the  probe  was  removed  and  transferred  to  an 
adjacent  location. 

In  each  experiment,  the  gas-sampling  probe  was  inserted  into  the  exhaust  gas 
stream  in  such  a way  to  assure  probing  the  conditions  at  the  selected  exit  plane,  about 
0.533  m from  the  end  plate  and  0.508  m from  the  vertical  axis  of  the  combustion 
chamber  (see  Fig.  1).  The  extracted  gas  samples  were  analyzed  with  the  FTIR 
spectrometer  and  the  spectral  data  were  recorded.  These  spectra  were  used  to  identify 
the  species  present  in  the  reactor  emissions.  Separately  generated  calibration  spectra  of 
each  identified  species  was  then  used  to  quantify  the  species  in  the  gas  sample. 

3.  Results  of  Experiments 

3. 1.  Operating  and  Boundary  Conditions 

The  fuel  flow  rate  into  the  reactor  was  maintained  at  3.0  ± 0.02  kg  h’^  and  Table  3 
summarizes  the  estimated  uncertainties  in  the  measurement.  The  methanol  flow  rate 
was  measured  using  a COX  turbine  meter.  The  rotation  frequency  of  the  turbine  meter 
was  calibrated  as  a function  of  methanol  flow  rate,  and  the  analog  voltage  signal  from 
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the  turbine  meter  was  read  into  the  computer  used  for  data  acquisition  and  calibrated  as 
a function  of  the  turbine  frequency.  Both  of  the  uncertainties  associated  with  these 
calibrations  are  presented  in  the  table.  Combining  the  calibration  uncertainties  with  the 
standard  deviation  obtained  from  repeated  observations  yields  a combined  standard 
uncertainty  of  0.0157  kg  h'^ 

The  combustion  air  was  delivered  through  the  swirl  burner  at  a rate  of  56.7  ± 1 .7  m^ 
h'^  The  flow  rate  was  measured  using  a 6.35  mm  i.d.  COX  sonic  nozzle  for  which  the 
manufacturer  reports  a 3 % imcertainty.  This  uncertainty  is  significantly  larger  than 
those  associated  with  the  calibration  of  the  pressure  transducer,  the  uncertainty  of  the 
pressure  gauge,  or  the  random  errors  determined  from  repeated  observations.  Table  4 
summarizes  the  various  components  of  the  combined  uncertainty,  which  is  equal  to 
1.73  m^h'^ 

Thermocouples  (K-type)  were  used  to  measure  the  wall  temperature  and  the  gas 
temperature  at  the  exit  of  the  reactor.  The  locations  and  mean  temperatures  are 
summarized  in  Tables  5 and  6.  The  temperature  of  the  exit  gas  was  measmed  at  the 
same  locations  where  the  species  concentration  data  were  obtained.  The  uncertainty 
budget  for  the  temperature  measurements  is  presented  in  Table  7,  and  the  combined 
uncertainty  is  3 °C  and  5.5  °C  for  the  wall  and  gas  temperatures,  respectively. 
Calibration  of  the  sonic  nozzle  and  thermocouples  is  planned  prior  to  distribution  of  the 
final  report  for  the  baseline  case. 

3. 2.  Fuel  Spray  Measurements 

Figure  4A  presents  an  example  of  droplet  size  data  obtained  at  one  location  in  the 
spray  (z  = 5 mm,  r = 2.76  mm).  Size  distributions  (probability  as  a frmction  of  droplet 
diameter)  from  thirteen  experiments  are  shown  in  the  figure.  The  data  show  good 
agreement  over  the  thirteen  experiments,  with  the  largest  variation  in  the  probability 
density  function  occurring  near  a droplet  diameter  of  15  |im.  The  data  from  Fig.  4A  are 
shown  in  Fig.  4B  with  a cubic  spline  interpolation,  which  represents  the  mean  size 
distribution.  Figure  4C  presents  the  mean  distribution  obtained  from  the  cubic  spline 
interpolation  with  pseudo-confidence  bands  that  represent  an  approximate  95  % 
confidence  interval.  The  pseudo-confidence  bands  were  calculated  as  ± 2sn  , where  s 
is  the  standard  deviation  of  the  raw  data  and  n is  the  number  of  samples. 

Experimental  data  were  collected  at  seven  elevations  downstream  of  the  nozzle  (z  = 
5,  15,  25,  35,  45,  55,  and  65  mm).  These  data  are  presented  in  the  following  manner. 
Droplet  size  distributions  corresponding  to  numerous  radial  locations  at  each 
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downstream  location  (z-coordinate)  are  presented  in  Figs.  5 - 11;  axial  velocity 
distributions  (Vz)  are  presented  in  the  same  fashion  in  Figs.  12  - 18;  and  radial  velocity 
distributions  (Vr)  are  presented  in  Figs.  19  - 25.  The  figures  present  the  cubic 
interpolation  through  the  distributions.  The  pseudo-confidence  (95  %)  bands  are  also 
shown  in  the  figures. 

The  data  shown  in  Figs.  5 - 25  are  provided  on  the  disk  accompanying  this  report. 
The  figures  provide  a quick  reference  for  selecting  the  data  files  of  interest.  The  ASCII 
text  files  contain  four  columns  of  data.  The  four  columns  are  the  droplet  diameter  (or 
velocity),  the  cubic  spline  interpolation,  the  lower  confidence  band,  and  the  upper 
confidence  band,  respectively.  The  file  names  correspond  to  the  type  of  distribution 
and  location  in  the  spray.  Files  that  begin  with  D,  Z,  and  R correspond  to  size,  axial 
velocity,  and  radial  velocity  distributions,  respectively.  The  number  immediately 
following  the  letter  indicates  the  axial  position  downstream  of  the  nozzle,  and  the 
following  number  indicates  the  radial  coordinate.  For  example,  the  file  D35_2159.spl 
contains  the  size  distribution  measured  at  z = 35  mm  and  r = 21.59  mm.  The  radial 
coordinates  are  also  indicated  in  Figs.  5-25. 

Droplet  Sauter  mean  diameters,  D32,  are  presented  in  Fig.  26  for  each  of  the  seven 
downstream  locations.  The  mean  axial  and  radial  velocity  componenets  of  the  droplets 
are  presented  in  Figs.  27  and  28,  respectively.  In  addition,  the  mean  components  of  the 
droplet  velocities  are  summarized  in  the  files  axialvel.dat  and  radvel.dat  contained  on 
the  accompanying  diskettes.  The  files  contain  twenty-one  columns  that  correspond  to 
the  radial  coordinate,  velocity,  and  uncertainty  for  each  of  the  seven  elevations  shown 
in  the  figures. 

Figure  29  presents  droplet  number  densities  measured  with  the  PDI  at  the 
downstream  locations  corresponding  to  data  shown  in  Figs.  5-25.  The  volume  flux  of 
fuel  droplets  corresponding  to  the  same  locations  is  shown  in  Fig.  30.  The  number 
densities  and  volume  fluxes  are  corrected  for  the  instrument  response  time  and  rejected 
signals,  which  is  discussed  in  detail  elsewhere  (Widmann  et  al.,  1999b).  The  horizontal 
error  bars  express  the  uncertainty  in  the  radial  coordinate,  and  the  vertical  error  bars 
correspond  to  twice  the  standard  error  of  the  mean  (2sri  ).  Note  that  systematic  errors 
associated  with  the  instrument  are  not  included  in  the  figures  at  this  time  but  should  be 
considered  when  comparing  the  data  to  numerical  results.  Systematic  errors  associated 
with  PDI  are  discussed  further  below.  The  accompanying  diskettes  contain  the  data 
presented  in  Figs.  29  and  30  in  ASCII  format,  and  also  contain  the  Sauter  mean 
diameters.  The  seven  files  are  named  spray_l.dat,  spray_2.dat,  etc.,  and  correspond  to 


9 


different  elevations  downstream  of  the  nozzle.  The  coluirm  format  of  the  data  files  is 
summarized  in  Table  8. 

The  uncertainties  associated  with  droplet  number  density  and  flux  measurements 
deserve  additional  attention  due  to  the  difficulty  in  their  quantification  and  their 
importance  in  satisfying  mass  conservation.  Zhu  et  ah  (1993)  have  discussed  the 
difficulties  of  measuring  number  densities  and  fluxes  using  PDI.  Simultaneously 
measuring  the  number  density  and  volume  flux  is  particularly  difficult  because  the 
number  density  is  dominated  by  the  many  small  droplets  while  the  volume  flux  is 
strongly  dependent  on  the  larger  droplets.  Therefore,  if  the  size  range  of  the  droplets  in 
the  spray  is  greater  than  the  dynamic  range  of  the  instrument  (35:1  for  the  counter 
based  processors,  as  used  in  this  study),  the  simultaneous  measurement  of  number 
density  and  volume  flux  is  very  challenging.  Research  into  statistical  methods  of 
accurately  accounting  for  the  droplets  beyond  the  measurable  size  range  of  the  PDI  is  in 
progress  and  may  lead  to  improved  volume  flux  measurements. 

The  accuracy  of  both  number  density  and  volume  flux  measurements  depends  on 
several  factors  (e.g.,  the  dynamic  range  of  the  instrument  and  accurately  determining 
the  probe  volume).  There  are  three  main  factors  that  contribute  to  uncertainties  in  the 
probe  volume.  The  first  is  the  variation  of  the  probe  diameter  with  droplet  size,  which 
results  from  the  Gaussian  intensity  profile  of  the  laser  beam.  The  probe  volume  can  be 
approximated  as  a cylinder,  where  the  laser  beams  determine  the  diameter  and  the 
length  is  determined  by  the  100  pm  aperture  in  the  receiving  optics.  The  effective 
diameter  of  the  laser  beam  depends  on  the  droplet  size  because  larger  droplets  scatter 
significantly  more  light  than  smaller  droplets  (van  de  Hulst,  1951).  Smaller  droplets 
must  pass  through  the  center  of  the  probe  volume  where  the  laser  light  intensity  is  high 
while  larger  droplets  can  scatter  sufficient  light  to  be  detected  when  passing  through  the 
wings  of  the  Gaussian  profile.  Therefore,  the  probe  volume  is  larger  for  the  larger 
droplets  than  for  the  smaller  ones. 

The  second  factor  contributing  to  the  xmcertainty  in  the  probe  volume  is  the  error 
associated  with  measuring  its  length.  The  accuracy  with  which  this  dimension  can  be 
measured  depends  upon  the  resolution  of  the  receiver  optical  system.  Bachalo  et  al. 
(1988)  estimated  this  uncertainty  to  be  15  % for  a 100  pm  slit;  however,  they  note  that 
larger  particles  will  produce  greater  uncertainty  in  this  dimension.  The  third  factor  that 
can  contribute  to  the  uncertainty  in  the  probe  volume  is  the  particle  trajectory.  It  is 
known  that  when  the  droplets  have  a component  of  velocity  in  the  direction  of  the  laser 
beams  (the  direction  that  is  not  measured),  the  probe  volume  determined  by  the  PDI 
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software  will  be  in  error.  Unfortunately,  the  magnitude  of  this  error  is  difficult  to 
determine.  For  the  spray  investigated  in  this  study,  it  has  been  reported  that  the 
droplets  have  a non-negligible  tangential  component  of  velocity  near  the  nozzle,  but 
this  component  is  significantly  reduced  further  downstream  (Presser  et  al.,  1993).  This 
may  result  in  errors  in  number  density  and  flux  measurements  close  to  the  nozzle,  for 
example,  z < 20  mm. 

The  difficulties  associated  with  quantifying  the  measurement  uncertainties  is 
demonstrated  in  Fig.  31,  where  the  reported  number  density  is  given  at  three  locations 
within  the  spray.  The  data  correspond  to  two  instrument  settings.  In  one  case,  the 
droplet  size  and  two  components  of  the  velocity  are  measured,  while  in  the  other  case, 
the  size  and  only  one  component  of  the  velocity  are  measured.  The  measurements 
made  with  only  one  component  of  velocity  yielded  number  densities  5-7  times  as  large 
as  those  made  with  two  components  of  velocities.  Note  that  there  is  no  significant 
difference  in  the  validation  fi'action  between  experimental  runs  with  and  without  the 
radial  velocity  measurement.  The  validation  fi'action  is  the  fraction  of  measured  signals 
from  which  the  processor  was  able  to  determine  the  size  and  velocity  of  the  droplet,  and 
in  all  of  the  experiments  presented  in  Fig.  31  the  validation  fraction  was  -0.5.  Also, 
the  area  of  the  sample  volume  reported  by  the  PDI  software  was  the  same  for  both 
settings;  however,  it  will  be  shown  below  that  this  may  be  incorrect.  Figures  32  - 34 
reveal  possible  explanations  for  the  discrepancy  between  the  number  density 
measurements.  Figure  32  shows  size  distributions  obtained  with  and  vdthout  the  radial 
velocity  measurement.  Turning  off  the  second  velocity  component  results  in  a 
measured  size  distribution  that  is  shifted  toward  smaller  droplets.  It  is  possible  that  the 
inclusion  of  these  smaller  droplets  in  the  size  distribution  is  responsible  for  the 
substantial  increase  in  the  number  density. 

Figure  33  presents  data  corresponding  to  the  number  of  fiinges  that  each  droplet 
crossed  during  an  experimental  run  in  which  the  size  and  two  components  of  velocity 
were  measured.  The  closed  symbols  correspond  to  measurements  in  the  radial  direction 
(green  laser  beams),  and  the  open  symbols  correspond  to  measurements  in  the  vertical 
direction  (blue  laser  beams).  Note  that  there  are  open  symbols  that  are  obscured  by  the 
closed  symbols.  It  is  clear  from  Fig.  33  that,  in  general,  more  blue  fringes  are  crossed 
during  the  measurement  than  green  fringes.  This  suggests  that  the  probe  area  may  be 
larger  for  the  measurements  in  the  axial  direction  than  the  radial  direction.  This  would 
explain  the  discrepancy  in  Fig.  31  if  the  probe  area  determined  from  the  blue  laser 
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beams  were  used  in  calculating  number  density  and  volume  flux  when  both  axial  and 
radial  velocity  components  are  measured. 

Figure  34  presents  the  calculated  probe  area  for  measurements  with  and  without  the 
radial  velocity.  The  probe  area  was  determined  as  a function  of  droplet  size  from  the 
raw  data  shown  in  Fig.  33.  The  probe  area,  Ap,  can  be  calculated  from 


Ap(dp)  = Dp(dp)  Wsiit  / sin(e),  ( 1 ) 

where  dp  is  the  droplet  diameter.  Dp  is  the  laser  beam  diameter,  Wsut  is  the  slit  width, 
and  0 is  the  scattering  angle.  For  these  experiments,  the  slit  width  and  scattering  angle 
were  100  pm  and  30°,  respectively.  The  laser  beam  diameter  is  computed  from 

Dp(dp)  = 8 [N„axH)  - (2) 

where,  5 is  the  fringe  spacing  , which  is  determined  by  the  wavelength  of  the  laser  light 
and  the  intersection  angle  of  the  laser  beams.  For  this  experimental  configuration,  the 
fringe  spacing  was  4.00  pm  and  3.95  pm  for  the  green  and  blue  laser  beams, 
respectively.  The  processor  requires  a minimum  number  of  fringe  crossings,  Nmin,  for  a 
signal  to  be  validated.  Also,  there  will  be  a maximum  number  of  fringe  crossings, 
Nmax,  in  each  size  class  due  to  the  finite  size  of  the  laser  beams.  The  PDI  software 
records  the  number  of  fringe  crossings  for  each  droplet  measured;  therefore,  Nmin  and 
Nmax  can  be  computed  for  each  size  class  and  Eq.  (2)  can  be  used  to  compute  Dp  as  a 
function  of  droplet  size.  Figure  34  reveals  that  the  probe  area  is  different  for  the  two 
channels,  which  suggests  that  the  discrepancy  shown  in  Fig.  31  may  be  due  in  part  to 
the  PDI  using  a probe  area  that  is  calculated  incorrectly  for  one  of  the  channels. 

To  evaluate  the  accuracy  of  the  volume  flux  measurements,  the  mass  flow  rate  of 
fuel  is  determined  at  each  axial  position  by  integrating  over  all  values  of  the  radial 
coordinate,  r.  The  results  normalized  by  the  inlet  flow  rate  are  presented  in  Fig.  35  with 
some  values  obtained  from  the  literature  (Sellens,  1990).  The  experiment  conducted  by 
Sellens  was  not  performed  under  identical  conditions  as  those  reported  here,  but  the 
data  demonstrate  the  difficulty  in  measuring  volume  fluxes  close  to  the  nozzle  using 
PDI.  Note  that  the  volumetric  flow  rates  obtained  by  integrating  the  volume  flux  are 
much  lower  than  one  would  expect  near  the  nozzle,  but  further  downstream  the  results 
are  reasonable,  ranging  from  60  % to  80  % of  the  inlet  fuel  flow  rate.  A large 
tangential  component  of  the  droplet  velocity  near  the  nozzle,  for  example,  would 
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introduce  errors  into  the  calculation  of  the  probe  volume  (as  discussed  above),  and 
would  contribute  to  the  low  values  of  the  integrated  mass  flux  that  are  observed  in  Fig. 
35.  Further  downstream,  where  the  tangential  component  of  the  velocity  has  decreased, 
the  error  would  not  be  as  pronounced.  However,  similar  results  were  reported  by 
Bulzan  et  al.  (1992)  for  non-swirling  conditions,  which  indicate  that  if  the  tangential 
component  of  velocity  is  contributing  to  the  error  it  is  not  the  only  important  factor. 

If  the  data  obtained  in  this  study  are  to  be  used  for  input  and  validation  of 
multiphase  combustion  models,  the  lower  values  of  the  droplet  volume  fluxes  shown  in 
Fig.  30  must  be  corrected  for  these  effects.  The  typical  manner  in  which  this  is 
accomplished  is  to  scale  the  number  densities  and  fluxes  presented  in  Figs.  29  and  30 
so  that  the  mass  flow  rate  at  z = 5 mm  is  consistent  with  the  inlet  fuel  flow  rate.  This 
procedure  appears  reasonable,  and  is  understandable  when  one  considers  the  lack  of 
alternatives.  However,  we  note  that  this  scaling  procedure  may  introduce  significant 
errors  into  the  assumed  spray  distribution.  If  this  procedure  is  adopted  to  specify  inlet 
conditions  for  CFD  models,  it  may  be  impossible  to  obtain  satisfactory  agreement 
between  the  computational  results  and  the  experimental  data.  Figure  36  illustrates  the 
difficulties  associated  with  this  procedure.  The  number  distribution  and  mass 
distribution  obtained  during  one  experimental  run  are  shown  in  the  figure.  While  the 
experimental  data  appear  to  capture  the  number  distribution  satisfactorily,  the  mass 
distribution  is  clearly  truncated  for  diameters  greater  than  60  |J.m.  The  polydispersity  of 
the  spray  makes  it  difficult  to  accurately  measure  both  number  and  mass  distributions, 
and  much  better  results  are  obtained  for  monodisperse  sprays.  Unfortunately,  an 
accurate  representation  for  the  number  distribution  that  accounts  for  all  of  the  fuel 
entering  the  reactor  is  critical  for  CFD  models.  Scaling  the  measured  fluxes  as  is 
commonly  done  may  result  in  overestimating  the  number  of  small  droplets  because  a 
few  large  droplets,  which  can  contain  a non-negligible  mass  fi-action  of  the  fuel,  are 
neglected. 

In  this  report,  spray  flame  measurements  are  presented  for  the  purpose  of  validating 
spray  combustion  models  and  submodels.  However,  it  is  necessary  to  quantify  the 
measurement  uncertainties  if  these  data  are  to  be  compared  with  multiphase  combustion 
models.  The  expressed  uncertainties  follow  the  guidelines  put  forth  by  NIST  (Taylor 
and  Kuyatt,  1994)  and  the  National  Conference  of  Standards  Laboratories  (1998). 
Uncertainties  are  classified  as  random  (Type  A)  or  systematic  (Type  B),  and  it  is 
assumed  that  the  combined  variance  (s  ) is  the  sum  of  the  component  variances. 
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There  are  many  factors  that  contribute  to  the  uncertainties  of  PDI  measurements 
such  as  uncertainties  associated  with  the  beam  separation,  transmitting  lens  focal 
length,  receiving  lens  focal  length,  scattering  angle,  and  others.  Sellens  (1990) 
summarized  the  various  contributions  and  provided  quantitative  estimates  for  his 
system.  Sellens  concluded  that  the  systematic  uncertainties  for  the  system  he  used  were 
4 % and  3 % for  the  velocity  and  size  measurements,  respectively.  These  uncertainties 
are  lower  than  those  reported  by  Bulzan  (1995),  who  used  a system  similar  to  the  one 
used  for  the  experiments  reported  here.  Bulzan  estimated  the  drop  size  uncertainty  in 
an  isothermal  monodisperse  drop  stream  to  be  6.5  %,  and  noted  that  it  would  probably 
be  larger  in  an  evaporating,  polydisperse  spray.  McDonell  and  Samuelsen  (1990) 
reported  that  the  size  error  due  to  variation  in  the  refractive  index  of  a methanol  spray 
when  the  temperature  was  varied  from  room  temperature  to  the  boiling  point  was 
3.5  %.  Because  the  experiments  described  here  were  conducted  with  methanol,  and 
with  the  same  scattering  angle  as  that  used  by  McDonell  and  Samuelsen,  we  will  use 
this  estimate  of  the  uncertainty.  Bulzan  reported  an  estimated  uncertainty  in  the 
velocity  measurements  of  10  %.  Taylor  et  al.  (1994)  examined  the  errors  associated 
with  PDI  measurements  using  monodisperse  polystyrene  latex  (PSL)  spheres  ranging  in 
diameter  from  0.7  |i.m  to  1 0 fim.  They  found  the  systematic  error  of  the  instrument  to 
be  a strong  function  of  particle  size,  with  errors  greater  than  50  % for  dp  ==  0.7  jiim  and 
errors  in  the  range  7 % - 14  % for  PSL  spheres  with  dp  > 1 jim. 

Droplet  number  density  and  flux  measurements  are  the  most  difficult  to  make 
because  they  require  accurate  measurements  of  droplet  size,  droplet  velocity,  and  probe 
volume.  Also,  the  number  density  and  volume  flux  measurements  are  strongly 
dependent  on  the  small  and  large  droplets,  respectively.  The  probe  volume  is  a 
function  of  the  droplet  diameter  due  to  the  Gaussian  intensity  profile  of  the  laser  beams 
as  discussed  above.  The  PDI  software  attempts  to  correct  for  the  size  dependence; 
however,  the  manufacturer  reports  an  uncertainty  of  7.5  % in  determining  the  sample 
volume  (Aerometrics,  1987).  In  addition,  multiple  particles  can  occupy  the  probe 
volume  simultaneously  which  "will  result  in  the  signal  being  rejected  and  the  number 
density  being  under-reported.  This  is  not  expected  to  be  a significant  source  of  error  in 
the  measurements  presented  in  this  report  because  the  number  of  rejections  was  very 
small  (<  0.5  %)  when  the  sizing  was  turned  off  and  only  the  velocity  components  were 
measured.  If  multiple  droplets  occupied  the  probe  volume  simultaneously  the  signal 
would  have  been  rejected  during  these  velocity  measurements. 
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When  using  diagnostic  techniques  based  upon  light  scattering  in  flames,  there  is  the 
possibility  that  local  variations  in  the  refractive  index  due  to  temperature  and 
concentration  variations  can  result  in  beam  steering.  The  PDI  manufacturer  has 
investigated  beam  steering  in  flames  and  reports  that  the  effect  on  the  measurements  is 
not  critical  (Taylor  et  ah,  1994).  Beam  steering  can  result  in  a loss  of  signal  so  that  the 
number  density  and  volume  flux  measurements  may  be  too  lo'w;  however,  this  is  not 
expected  to  introduce  significant  errors  into  the  size  and  velocity  measurements. 

Table  9 summarizes  the  uncertainty  budget  for  the  PDI  measurements.  There  are 
several  points  that  deserve  attention.  First,  the  random  errors  associated  with  repeated 
measurements  (Type  A)  are  much  less  than  the  systematic  errors  (Type  B).  The 
repeatability  of  the  measurements  reported  in  this  paper  was  very  good,  which  resulted 
in  consistent  results  during  the  several  month  period  while  data  were  collected. 
Unfortunately,  the  systematic  errors  associated  with  the  phase  Doppler  interferometer 
are  much  larger  than  the  random  errors.  The  result  is  that  although  the  spray  was 
consistent  from  one  day  to  another,  the  ability  to  quantify  the  droplet  size,  velocity, 
number  density,  etc.  is  limited  by  the  instrument.  Second,  it  should  be  noted  that  the 
uncertainties  associated  with  the  number  density  and  volume  flux  measurements  are 
much  greater  than  those  associated  vvdth  the  size  and  velocity  measurements.  This  is 
true  of  both  Type  A and  Type  B errors. 

3.S.  Exhaust  Chemical  Species  Measurements 

The  absorption  spectra  obtained  with  the  FTIR  spectrometer  indicate  that 
combustion  was  incomplete.  Major  species  identified  by  the  FTIR  data  include  CO2, 
CO  and  CH3OH.  No  minor  components  or  reaction  intermediates  were  identified 
perhaps  due  to  the  short  path  length.  Figure  37  shows  the  mole  fraction  of  CO2,  CO, 
and  CH3OH  at  the  exit  plane  of  the  reactor  (see  Fig.  1).  The  species  concentrations  are 
fairly  uniform  across  the  reactor  exit  plane,  suggesting  that  there  is  good  mixing  in  this 
region.  Also,  the  concentration  of  CO2  is  approximately  50  times  greater  than  that  of 
CO,  which  indicates  that  the  rate  of  methanol  conversion  to  CO2  and  H2O  is  occurring 
at  a rate  approximately  50  times  faster  than  the  conversion  to  CO  and  H2O.  The 
uncertainty  budget  for  the  species  concentration  measurements  is  presented  in  Table  1 0. 
The  concentrations  of  NO2  and  N2O  were  too  low  to  detect  with  the  FTIR 
measurements.  A NOx  analyzer  is  being  installed  to  measure  the  concentrations  of 
these  chemical  species. 
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Figure  38  presents  theoretical  species  concentrations  within  the  reactor  emissions. 
The  calculations  were  performed  assuming  that  only  two  reactions  were  significant. 
These  are  the  conversion  of  CH3OH  to  CO2  and  the  conversion  of  CH3OH  to  CO.  The 
rate  of  CO2  production  is  assumed  to  be  50  times  faster  than  the  production  of  CO, 
which  was  shown  above  to  be  consistent  with  the  FTIR  data.  The  filled  circles  in  Fig. 
38  show  the  measured  concentrations  of  CO,  CO2,  and  CH3OH.  These  are  spatially 
averaged  values  obtained  from  Fig.  37.  The  velocity  profiles  were  not  required  for 
these  calculations  because  of  the  spatially  uniform  species  concentrations.  The 
comparison  of  experimental  data  and  theoretical  calculations  in  Fig.  38  show  that  the 
conversion  of  CH3OH  is  ~ 80  %,  and  that  mass  is  conserved  in  the  system. 

4.  Summary 

Experimental  data  with  uncertainty  budgets  were  obtained  from  a reacting 
methanol  spray  under  well-characterized  conditions  for  the  purpose  of  validating 
multiphase  combustion  models  and  submodels.  Input  and  boundary  conditions  are 
provided  for  a defined  baseline  case.  Droplet  size  and  velocity  distributions,  number 
densities,  and  volume  fluxes  were  measured;  however,  the  uncertainty  associated  with 
the  number  density  and  volume  flux  measurements  are  considered  high.  Gas-phase 
species  measurements  were  made  using  FTIR  spectroscopy,  and  provide  exhaust 
concentrations  of  CO,  CO2,  and  CH3OH.  The  species  measurements  indicate  that  the 
conversion  of  methanol  in  the  reactor  is  ~ 80  %.  Wall  temperature  data  are  provided  as 
a function  of  axial  position,  as  well  as  the  temperature  of  the  exhaust  gas.  Gas-phase 
velocity,  temperature,  and  heat  flux  measurements  are  planned  for  the  next  stage  of  this 
investigation. 
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Table  1.  Operating  conditions  for  the  baseline  case. 


Fuel  Type 

Methanol 

Fuel  Flow  Rate 

3.0  kgh'^ 

Fuel  Temperature 

Ambient 

Equivalence  Ratio 

0.3 

Air  Flow  Rate 

56.7  m^  h’^ 

Air  Temperature 

Ambient 

Vane  Angle 

o 

o 

Swirl  Number 

0.47  (geometrical  correlation) 

0.54  (CFD  results) 

Flame  Standoff  Distance 

~ 5 mm 

Chamber  Pressure 

Ambient 
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Table  2.  Uncertainty  budget  for  the  location  of  PDI  probe  volume. 


Source  of 
uncertainty 

Standard 

uncertainty 

(mm) 

Type  B 

axial  position,  Az 

0.17 

horizontal  position.  Ax 

0.083 

horizontal  position.  Ay 

0.17 

spray  angle,  a6 

0.33® 

radial  coordinate,  Ar 

^ /•  + z tan( A0)  j 
/•-ztan(A0) 
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Table  3.  Uncertainty  budget  for  the  fuel  flow  rate  (3.0  kg  hr'*). 


Source  of 
uncertainty 

Standard 

uncertainty 

(kghr-‘) 

Type  A 

Type  B 

fuel  flow  / turbine  frequency 
calibration 

0.0042 

turbine  meter  / DAQ  signal 

0.0020 

repeated  observations 

0.015 

Combined  standard  imcertainty:  0.0157  kg  hr'* 

22 


Table  4.  Uncertainty  budget  for  the  combustion  air  flow  rate  (56.7  m hr'  ). 


Source  of 
uncertainty 

Standard 
uncertainty 
(m^  hr'*) 

Type  A 

T5T)e  B 

flow  rate  / pressure  calibration 
reported  by  manufacturer 

1.70 

pressure  / voltage  calibration 

0.252 

comparison  with  calibrated 
pressure  gauge 

0.170 

repeated  observations 

0.024 

Combined  standard  uncertainty:  1 .73  m^  hr'* 
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Table  5.  Wall  temperatures. 


Thermocouple 

Height 

(mm) 

Mean 

Temperature 

(°C) 

TC#1 

99 

93 

TC  #2 

165 

93 

TC  #3 

231 

95 

TC#4 

297 

99 

TC  #5 

429 

103 

TC  #6 

561 

112 

TC  #7 

693 

114 

TC  #8 

826 

121 

TC  #9 

958 

133 

TC#10 

1090 

136 
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Table  6.  Exit  gas  temperatures 


Radial 

Coordinate,  p 
(mm) 

Polar 
Angle,  0 
(deg) 

Mean 

Temperature 

ro 

0 

- 

244 

44.5 

0 

265 

88.9 

0 

282 

133.4 

0 

288 

177.8 

0 

266 

44.5 

90 

251 

88.9 

90 

258 

133.4 

90 

271 

177.8 

90 

277 

44.5 

180 

221 

88.9 

180 

207 

133.4 

180 

182 

177.8 

180 

160 
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Table  7.  Uncertainty  budget  for  the  temperature  measurements. 


Source  of 
uncertainty 

Standard 

imcertainty 

(°C) 

Type  A 

Type  B 

uncertainty  reported 

by  the  manufacturer 

2.2 

repeated  observations 

wall  temperature 

<2.0 

gas  temperature 

<5.0 

Combined  standard  uncertainty: 

wall  temperature 

3.0 

gas  temperature 

5.5 
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Table  8.  ASCII  data  file  format 


File  Name 

Axial  Coordinate,  z 
(mm) 

spray_l.dat 

5 

spray_2.dat 

15 

spray_3.dat 

25 

spray_4.dat 

35 

spray_5.dat 

45 

spray_6.dat 

55 

spray_7.dat 

65 

Column 

Variable 

1 

Axial  Coordinate,  z (mm) 

2 

Radial  Coordinate,  r (mm) 

3 

Ar  (mm) 

4 

Sauter  Mean  Diam,  D32  (jim) 

5 

AD32  (urn) 

6 

Number  Density,  Nd  (cm*^) 

7 

ANd  (cm'^) 

8 

Volume  Flux,  Fv  (cm^  cm'^  s'^) 

9 

AFv  (cm^  cm*^  s'^) 

10 

Spray  Angle,  0 (degrees) 
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Table  9.  Uncertainty  budget  for  the  PDI  measurements. 


Source  of 
uncertainty 

Standard 

uncertainty 

(%) 

Type  A 

Type  B 

droplet  size  measurement: 

repeated  observations 

1.0 

uncertainty  reported  in  literature 

7.4 

droplet  velocity  measurement: 

repeated  observations 

1.5 

imcertainty  reported  in  literature 

10 

number  density  measurement: 

repeated  observations 

9.1 

volume  flux  measurement: 

repeated  observations 

6.7 

Combined  standard  uncertainty: 
droplet  size  7.45 

droplet  velocity  10.1 
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Table  10.  Uncertainty  budget  for  the  species  concentration  measurements. 


Source  of 
uncertainty 

Standard 

uncertainty 

(ppm) 

Type  A 

Type  B 

calibrations 

CH3OH 

5 X 10'^ 

CO 

10-^ 

CO2 

10'^ 

repeated  observations 

CH3OH 

4x  10^ 

CO 

7x  10'^ 

CO2 

4x  10'^ 

Combined  standard  uncertainty 


CH3OH 

4x  10"^ 

CO 

1.2  X 10"* 

C02 

4.1  X 10'^ 
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Figure  1.  Schematic  of  (A)  the  reference  spray  combustion 
facility  and  (B)  the  nozzle  used  to  generate  the  fuel 
spray. 
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Figure  2.  Relevant  dimensions  of  the  reactor  and  the 

burner  used  in  the  experiments  described  here, 
(dimensions  in  millimeters) 


Figure  3.  Schematic  of  the  image  observed  through  the  viewing 
port  of  the  PDI  receiving  optics  when  the  probe  volume 
is  located  at  (A)  x = y = 0 and  z = 1 mm,  and  (B)  x = y = 0 
and  z = 3 mm. 
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PROBABILITY 


Figure  4.  Droplet  size  distributions  at  z = 5 mm  and  r = 2.76  mm.  Presented 
are  (A)  the  raw  data,  (B)  the  raw  data  with  the  cubic  spline,  and 
(C)  the  cubic  spline  with  confidence  bands. 
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Figure  12.  Droplet  axial  velocity  distributions  at  various  radial 
coordinates  5mm  downstream  of  nozzle 
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Figure  13.  Droplet  axial  velocity  distributions  at  various  radial 
coordinates  15mm  downstream  of  nozzle 
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Figure  14.  Droplet  axial  velocity  distributions  at  various  radial 
coordinates  25mm  downstream  of  nozzle 
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AXIAL  VELOCITY,  m/s 


Figure  15.  Droplet  axial  velocity  distributions  at  various  radial 
coordinates  35mm  downstream  of  nozzle 


Ainiavaoud 

44 


AXIAL  VELOCITY,  m/s 


Figure  16.  Droplet  axial  velocity  distributions  at  various  radial 
coordinates  45mm  downstream  of  nozzle 
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AXIAL  VELOCITY,  m/s 


Figure  17.  Droplet  axial  velocity  distributions  at  various  radial 
coordinates  55mm  downstream  of  nozzle 
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AXIAL  VELOCITY,  m/s 


Figure  18.  Droplet  axial  velocity  distributions  at  various  radial 
coordinates  65mm  downstream  of  nozzle 
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AXIAL  VELOCITY,  m/s 


Figure  19.  Droplet  radial  velocity  distributions  at  various  radial 
coordinates  5mm  downstream  of  nozzle 
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Figure  20.  Droplet  radial  velocity  distributions  at  various  radial 
coordinates  15mm  downstream  of  nozzle 
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RADIAL  VELOCITY,  m/s 


Figure  21.  Droplet  radial  velocity  distributions  at  various  radial 
coordinates  25mm  downstream  of  nozzle 
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RADIAL  VELOCITY,  m/s 


Figure  22.  Droplet  radial  velocity  distributions  at  various  radial 
coordinates  35mm  downstream  of  nozzle 
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RADIAL  VELOCITY,  m/s 


Figure  23.  Droplet  radial  velocity  distributions  at  various  radial 
coordinates  45mm  downstream  of  nozzle 
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RADIAL  VELOCITY,  m/s 


Figure  24.  Droplet  radial  velocity  distributions  at  various  radial 
coordinates  55mm  downstream  of  nozzle 
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Figure  25.  Droplet  radial  velocity  distributions  at  various  radial 
coordinates  65mm  downstream  of  nozzle 
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SAUTER  MEAN  DIAMETER,  urn  SAUTER  MEAN  DIAMETER, 


Figure  26.  Sauter  mean  diameter,  D32,  at  seven  locations  downstream 
of  the  nozzle  (z  = 5,  15,  25,  35,  45,  55,  65  mm). 
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Figure  27.  Droplet  axial  velocity,  v^,  at  seven  locations  downstream 
of  the  nozzle  (z  = 5, 15,  25,  35,  45,  55,  65  mm). 
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Figure  28.  Droplet  radial  velocity,  v^,  at  seven  locations  downstream 
of  the  nozzle  (z  = 5, 15,  25,  35,  45,  55,  65  mm). 
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Figure  29.  Droplet  number  densities  at  seven  locations  downstream 
of  the  nozzle  (z  = 5,  15,  25,  35, 45,  55,  65  mm). 
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VOLUME  FLUX,  cm  cm'  s'  VOLUME  FLUX,  cm^  cm  ' 


Figure  30.  Droplet  volume  fluxes  at  seven  locations  downstream 
of  the  nozzle  (z  = 5,  15,  25,  35, 45,  55,  65  mm). 
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NUMBER  DENSITY,  cm' 


Figure  31.  A comparison  of  reported  number  densities  at  three  locations 

in  the  spray  with  one  and  two  components  of  velocity  measured. 
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Figure  32.  Comparison  of  reported  size  distributions  with  and  without 

the  radial  component  of  velocity  measured  at  (A)  z = 5 mm,  and 
(B)  z = 25  mm. 
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FRINGE  COUNT 


Figure  33.  Raw  data  obtained  with  the  PDI  showing  the  number  of  fringes 
crossed  for  each  measurement  direction  (axial  and  radial). 
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PROBE  AREA, 


Figure  34.  The  calculated  probe  area  based  upon  the  number  of  Singes  that 
individual  droplets  crossed. 
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NORMALIZED  MASS  FLOW 


Figure  35.  The  normalized  mass  flow  rate  of  fuel  at  various  downstream 

locations.  Values  from  the  literature  are  shown  for  comparison. 
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Figure  36.  Number  and  mass  distributions  measured  with  the  PDI  illustrating 
the  difficulty  of  measuring  the  number  and  mass  distributions 
simultaneously. 
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MOLE  FRACTION 


Figure  37.  The  exit  concentrations  of  various  species  in  the  reactor. 
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MOLE  FRACTION 


Figure  38.  Theoretical  species  concentrations  in  the  reactor  emissions. 
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